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Abstract 

We investigate the effect of dark energy on the limits on the shear 
anisotropy in spatially homogeneous Bianchi cosmological models ob- 
tained from measurements of the temperature anisotropics in the cos- 
mic microwave background. We shall primarily assume that the dark 
energy is modelled by a cosmological constant. In general, we hnd that 
there are tighter bounds on the shear than in models with no cosmolog- 
ical constant, although the limits are (Bianchi) model dependent. In 
addition, there are special spatially homogeneous cosmological models 
whose rate of expansion is highly anisotropic, but whose cosmic mi- 
crowave background temperature is measured to be exactly isotropic 
at one instant of time. 
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1 Introduction 



Observations suggest that the Universe has 0.98 < Vt t ot < 1-08 with radiation 
(R), baryons (B), dark matter (DM) and dark energy (DE) each contributing 
n fi ~5x 10- 5 ,fi B ~ 0.04, fi DM ~ 0.26,fi DjB ~ 0.7, respectively [1], where 
the variables Qi = Pi/p c give the fractional contribution of different compo- 
nents of the Universe (i denoting baryons, dark matter, radiation, etc.) to the 
critical energy density p c = 3Hq/8tiG, where H = (a/a) is the present Hub- 
ble expansion rate. Consequently, we have that f2 = flu + Qb + &dm ~ 0.3. 
The simplest model for a fluid with negative pressure is the cosmological con- 
stant [2] whence nearly seventy per cent of the energy density in the Universe 
is unclustered and exerts negative pressure; f^o = &de ~ 0.7. 

It is of interest to study the effect of dark energy on the cosmological limits 
from cosmic microwave background (CMB) measurements on the shear (and 
the rotation, etc.), especially in spatially homogeneous (SH) cosmological 
models. In the presence of the cosmological constant, an SH universe that is 
very nearly flat now and that was very nearly flat at the time of last scattering 
can become significantly negatively curved at intermediate periods. At earlier 
times the matter dominates and the curvature is dynamically negligible, while 
at later times fix can dominate. For a currently almost flat universe with 
parameters close to the observed values of fl\ ~ 0.7 and Q ~ 0.3, the 
curvature is at most only a few percent. At intermediate times, however, the 
curvature can be dynamically significant. 

The observed temperature anisotropy (from COBE and WMAP [3]) of 
the CMB radiation, AT/T, on large angular scales is determined by the value 
of the shear at the surface of last scattering of the radiation. To evaluate 
large-scale anisotropies we must evaluate the peculiar redshift a photon will 
feel from the epoch of last scattering (Is) until now (0) by integration along 
null geodesies in the background spacetime. 

The usual paradigm invoked to understand cosmological observations, 
motivated by inflation and consistent with the smoothness of the CMB, is 
that the Universe was spatially homogeneous and isotropic spacetime at early 
times. The theory is then that some physical mechanism (such as inflation) 
generated perturbations that subsequently evolved through gravitational col- 
lapse to form the structures we now observe. Within this paradigm, the 
Universe is closely modelled by a Friedmann-Robertson- Walker (FRW) cos- 
mology at all times since the earliest epoch (and particularly since recombi- 
nation). In such models it is often the case that the integrated effect from a 
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small shear is not important. 

We utilize dynamical systems methods using expansion-normalized vari- 
ables in SH cosmological models [I], generalized to include a cosmological 
constant, to revisit the cosmological constraints on the Hubble- normalized 
shear scalar X = a/(^/3H) (for example) from CMB anisotropy measure- 
ments. 

2 Background 

Bianchi models provide a description of general SH anisotropic cosmologies. 
There are distinct features depending on the overall geometry and homo- 
geneity class of the model |5|. In a pioneering paper, Collins and Hawking 
used analytical arguments to find upper bounds on the amount of shear (and 
vorticity) in the Universe today, from the absence of any detected CMB 
anisotropy [6]. A detailed numerical analysis of such models [7] used experi- 
mental limits on the dipole and quadrupole to refine limits on the universal 
shear (and rotation). More recently, it has been argued [8] that there is no 
'isotropy problem' in Bianchi type VII^ cosmological models in which shear 
(and vorticity) have decayed only logarithmically since the Planck time, in 
that the present amplitude of CMB fluctuations is compatible with current 
observational limits in these models. 

The limits could be strengthened in a number of ways; for example, if we 
assume particular classes of cosmological models (i.e., spatially homogeneous 
models) for which exact evolution laws are known we can replace assumptions 
and approximations by exact relations. For particular Bianchi models with 
specific evolution equations, stronger limits are indeed possible. In particular, 
in general shear anisotropy grows with time (i.e., Bianchi models do not 
generically isotropize). However, for very special models (i.e., not generic) of 
Bianchi type I and V the anisotropy decays, and hence much stronger limits 
on present day shear may be obtained in these models. 

In particular, in [9] it was shown that small quadrupole anisotropics in 
the CMB imply severe limits on spacetime anisotropy in a Bianchi type I 
model (assumed to be a small perturbation of a FRW model). Using an 
exact Bianchi I cosmological solution, the quadrupole component of CMB 
temperature distribution found by COBE |3J (expressed as an upper limit 
on AT/T), implies that So < 2 x 10~ 9 . Similarly, strong limits can be ob- 
tained in Bianchi type V models (and a special class of spatially homogeneous 
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Bianchi type Vllh models [8]). However, such severe limits are not generic. 

The above limits on the shear (and rotation) are not based on 'full sky' 
maps; only the simple quadrupole (and not higher multipoles, which produce 
hotspots/ spirals) is utilized. In general SH models the distorted quadrupole is 
combined with spiral geodesic motion, and by solving the geodesic equations 
describing the motion of photons numerically, limits can be derived from 
the whole COBE sky map. We note that unlike the open type VII^ models 
(for example), the closed Bianchi type IX models exhibit neither geodesic 
focusing nor the spiral pattern, even in the presence of vorticity [7]. 

A special class of Bianchi type Vllh models, which have a logarithmic 
dynamical evolution and E that is effectively 'frozen in' at late times, were 
studied in [7j and it was concluded that limits on the shear are as strong 
as | So | ~ 1CT 9 in the real Universe. In further work [10], the large-scale 
cosmic microwave background anisotropies in spatially homogeneous, glob- 
ally anisotropic cosmologies were investigated, and improvements on previous 
bounds on the total shear in the Universe were obtained by performing a sta- 
tistical analysis to constrain the allowed parameters of a Bianchi model of 
type Vllh using the cumulative data from COBE. Consequently, very strong 
upper limits on the amount of shear (and vorticity) were obtained; limits 
which are typically one to two orders of magnitude higher than constraints 
relying entirely on the quadrupole. Moreover, in discarding information from 
higher moments, the comparison is not sensitive to the small-scale structure 
present in anisotropic models that is associated either with the spiral pattern 
or with geometrical focusing when Q < 1. This bound ([10]) can be improved 
upon by a factor of \/3 by using different statistics [TT] . 

If the CMB fluctuations about isotropy are small, then so are deviations 
from spatial homogeneity and isotropy [T2], provided that the Universe after 
last scattering can be modelled as exactly dust (plus collision-free radia- 
tion, and may include a cosmological constant). Model-independent limits 
on the present day strengths of large-scale shear in the Universe in the im- 
portant case of an almost isotropic CMB, which also take into account inho- 
mogeneities as well (assuming constraints on the size of temporal and spatial 
derivatives), can then be obtained. A general limit on the present value of 
S can be obtained using the quadrupole and octopole; S < 10~ 4 [12]. This 
general limit is a much weaker limit on the shear than discussed above, which 
implies that the Bianchi VII^ models (for example) are too special to draw 
conclusions from because of their atypical evolution. However, it was shown 
[T2] that stronger limits on the shear [7J can be recovered in a (Bianchi VILJ 
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model which is Bianchi I to first order |12j . 

It is worth mentioning that since the Einstein equations are highly non- 
linear, the limit as the (anisotropic) curvatures tend to zero (in some ap- 
propriate sense) in Bianchi models does not necessarily produce the same 
results as those in a Bianchi I model with zero curvature. Although limits 
can often be strengthened by considering a particular class of Bianchi cos- 
mological models for which exact evolution equations are known, the results 
are not generic. Indeed, for very specific models of Bianchi type I and V 
which isotropize, very severe limits are obtained. In general, specific limits 
are obtained in each class of SH models, which include integrated and non- 
integrated effects, which combine distorted quadrupole and spiral geodesic 
motion, and which include limits from numerical and statistical analyses de- 
rived from the full COBE sky map. 



3 Analysis 

We shall investigate anisotropic cosmological models with pressureless matter 
(dust) and a cosmological constant (assumptions which are justified since 
recombination). Typically, if the shear E is small throughout the evolution 
since recombination, the integrated effect on the CMB is not important. 
However, if the shear is not small the effect can be important. 

In the case of a Bianchi type I model, and exact solution (for dust and a 
cosmological constant) is again possible. We find that 

V 2.-3* y 2„-3t 



s 2 e- 3 * + n Ao e 3t + n n Ao e 3t + n ' 

where So (for example) is the value of the shear at the present logarithmic 
time t — 0, and the constants satisfy 

n + n Ao + Eg = i. (2) 

The approximation can be integrated exactly to give 
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where we will use t[ s = —7 as the logarithmic time of last scattering!]] Im- 
posing the quadrupole limitj^] 

r° 

/ Y,dt < 1(T 5 , (4) 

we then obtain 

S < 4.13 x KT^v 7 ^. (5) 

Therefore, the limit is strengthened (for a smaller value of Qq, corresponding 
to a larger Qa ) over the case with no cosmological constant (for example, by 
a factor of about 2 if Qq = 0.3). The physical reason for this strengthening 
is that for a more dominant cosmological constant (as measured by Q\ ), 
the Hubble-normalized shear E tends to zero faster (see equation ([1])). As a 
result, given a fixed amount of integrated shear (e.g., 10~ 5 ), the present day 
value So is smaller in the presence of a dominant fl\ . 

The limit on the shear in Bianchi I models from the quadrupole is strength- 
ened when a non-zero cosmological constant is present. There will be similar 
results for Bianchi type V models, which also isotropize to the future (and, in 
addition, for Bianchi type VII^ models). However, for general models which 
do not isotropize to the future (without the presence of dark energy), the 
limits on the shear will become much much tighter. 

1 The logarithmic time t is defined as the logarithm of the length scale: t = ^o et_tt \ 
where we have set to — above. The logarithmic time is therefore related to the cosmic 
time through dt dt — H, where H is the Hubble scalar. To determine the value of t/ s , we 
take the temperature of the last scattering surface to be 3000K. Using the linear relation 
between the CMB temperature and i -1 , we obtain 

t la = In f^j = - In f^j = - ln(3000/2.7) w -7. 

The value of ti s can be adjusted for different values of the temperature of the last scattering 
surface. 

2 As an idealization, we assume that the last scattering temperature is isotropic and 
that all the quadrupole anisotropy in the CMB comes from the integrated shear. We use 
the upper bound ^ as a rough (order of magnitude) approximation of the bound on the 
quadrupole, a2 < 10~ 5 , in order to forgo numerically calculating the CMB map. The 
integral in (TJ| is equal to 

I ^G^cosmic 

•Atcosmic) ls V3 

when expressed as an integral of the shear <j over cosmic time. 
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As an illustration, let us consider the locally-rotationally symmetric (LRS) 
Bianchi type III cosmological models. The evolution equations for expansion- 
normalized variables can be obtained by setting k — 1/2 in [T3| Section 2.1.2] 
and including a cosmological constant: 

E' x = (g-2)S x -fi fe (6) 

fi' fc = 2(g + £ x )fi fe (7) 

fi A = 2(g + l)fi A , (8) 

where 

q = 2£ 2 X + |fi - fi A (9) 

fi = 1 - S 2 X - fi fc - fi A . (10) 

Here, S x is the shear, fi& is the spatial curvature, fi A is the density parameter 
for the cosmological constant, fi is the density parameter for dust, and q is 
the deceleration parameter. 

Shown in Figured] is the state space of LRS Bianchi type III cosmological 
models with dust and a cosmological constant. The state space shown is a 
quarter of a unit sphere, with equilibrium points representing the flat FRW 
solution (F), the de Sitter solution (dS), the two Kasner solutions (T and Q) 
and the plane wave solution (P) (located at S x = — |). The FRW solution 
with a cosmological constant is represented by the vertical line F-dS. The 
state space is divided into regions A and B (separated by the surface F- 
P-dS), according to whether the shear S x changes sign over the course of 
evolution. 

The observed CMB temperature T is given in terms of an integral of the 
quantity s, which in this case only depends on E x and one direction cosine 
K 2 : 



T = Tis exp 



'1 + s) dt 



Is 

2\ 



S = -(l-3^)S x , (12) 
K' 2 = -3(1 - K\)K 2 ^. (13) 



Following the calculation in [TJ], the observed CMB temperature T is then 
given by 



e 2/3 



T = Ti e~ tls ^^^^^=^^^^^= (14) 

y/e** + (1 -e 6 A>)AT 2 2 (o)' 
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Figure 1: State space of LRS Bianchi type III cosmological models with dust 
and a cosmological constant. Arrows indicate the direction of evolution. The 
shaded surface indicates a set of initial conditions at present time yielding 
exactly isotropic CMB. 

where (3q is the integrated shear: 



If (3 = 0, then the CMB temperature is exactly isotropic. 

If the shear variable E x does not change sign (e.g., for initial conditions 
in the region B), then the viable models stay close to FRW during the whole 
evolution. Numerically integrating the null geodesic equations shows that, 
as expected, the limits on the shear are much tighter in the presence of a 
cosmological constant, tending to the Bianchi I constraints in the limit of 
negligible curvature. 

It is known that there are spatially homogeneous cosmological models 
whose CMB temperature is measured to be exactly isotropic (by all funda- 
mental observers), and hence indistinguishable from FRW, at one instant of 
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time, but whose rate of expansion is highly anisotropic [2] (recall that we 
only observe the CMB at one instant of time). These results do not contra- 
dict bounds on the shear found by [SI C3 E] , who restrict attention ab initio 
on spatially homogeneous models close to FRW. In particular, in the Bianchi 
type III models under investigation but with dust only (and no cosmological 
constant), if S x can change sign (i.e., for initial conditions in the region A 
[very close to the orbit F-P]), then the models need not stay close to FRW 
during the whole evolution and it is possible for the CMB to be exactly 
isotropic at one instant of time. From Figure HJ the curve at the base of the 
shaded surface shows that So can be very large (close to one half), while 
giving an exactly isotropic CMB. Initial conditions in a small neighbourhood 
of this curve gives a close-to-isotropic CMB. 

Similarly, in Bianchi type III models with dust and a positive cosmological 
constant, there is a surface of initial conditions yielding an exactly isotropic 
CMB and a small neighbourhood of this surface gives a close-to-isotropic 
CMB. The value of Sg on this surface is smaller than the value of Sq on 
the curve (i.e., adding a cosmological constant tightens the bound on S ). 
Note that the sign change in S x occurs (during the evolution) for all initial 
conditions on the indicated surface. As a result, imposing a limit on the 
CMB quadrupole will yield a looser bound on So than in Bianchi I models 
(or any model in which the shear components cannot change sign; a sign 
change can lead to a cancellation in the integral J Tidt). 

The same argument applies to non-LRS models, and generally in Bianchi 
models of all types, where there are sign changes in the shear [H] . However, 
complications in the analysis can arise because there may be no exactly 
isotropic surface to base the calculation on. In addition, obtaining constraints 
based upon the octopole may also pose problems, since J Hdt may be a poor 
approximation for the octopole. 

4 Discussion 

The WMAP data provide some of the most accurate measurements yet of 
the CMB and contribute to high accuracy determinations of cosmological 
parameters [3]. However, there are several studies that show that at large 
scales, the CMB is not statistically isotropic and Gaussian [15]. In [IS], a 
possible detection of a correlation between several independently discovered 
anomalies in the CMB sky measured by the WMAP data and a Bianchi Type 
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Vllh template [7J was reported. The predictions for the CMB anisotropy 
patterns arising in Bianchi type VII^ universes which include a dark energy 
component were presented in [17], and it was found that including a term 
Q\ > 0, can lead to the same observed structure as in the best-fit Bianchi 
type Vllh model of [16]. Although the best-fit Bianchi type Vllh model is 
not compatible with measured cosmological parameters, it was shown that 
removing a Bianchi component from the WMAP initial data release can 
account for several large angular scale anomalies and yield a corrected sky 
that is statistically isotropic. It was concluded that in the absence of an 
unknown systematic effect which could explain both the anomalies and the 
correlation, the WMAP data require an addition to the standard cosmological 
model that resembles the Bianchi morphology [T6] . 

We have only discussed orthogonal spatially homogeneous models here. 
Expanding SH perfect fluid cosmological models with a constant equation of 
state parameter 7 with a non-geodesic 'tilting' fluid congruence were stud- 
ied in |18j . It was shown that for stiff equations of state the tilt can give 
rise to an effective quintessential, or even phantom-like, equation of state. 
More importantly, it is of interest to investigate the effect of tilt on CMB 
radiation observations from the perspective of the observers moving with the 
fluid matter congruence and, in particular, to determine whether the tilt can 
significantly affect the constraints on the shear in SH models. 
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